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We study the persistence in a class of continuous stochastic processes that are stationary only 
under integer shifts of time. We show that under certain conditions, the persistence of such a 
continuous process reduces to the persistence of a corresponding discrete sequence obtained from 
the measurement of the process only at integer times. We then construct a specific sequence for 
which the persistence can be computed even though the sequence is non-Markovian. We show that 
this may be considered as a limiting case of persistence in the diffusion process on a hierarchical 
lattice. 
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I. INTRODUCTION 



In recent years, there has been a lot of interest in the 
study of persistence of fluctuations in different physical 
systems 0,0. Persistence P(t) is simply the probability 
that the deviation of the value of a fluctuating field from 
its mean value does not change sign up to time t. Persis- 
tence has been studied in many nonequilibrium systems 
|ij and also in diverse fields ranging from ecology || to 
seismology [|] . Theoretical studies include various mod- 
els of phase ordering kinetics M , diffusion equation HQ , 
reaction diffusion systems in pure Q as well as disordered 
environments fluctuating interfaces JTo| |, and various 
theoretical models |ll[] . Persistence or first passage prop- 
erties find simple applications in various chemical jL2|, 
biological [|l3| and granular systems In laboratory 

experiments, persistence has been measured in various 
experimental systems including breath figures |]l5f , liq- 
uid crystals ijlql , soap bubbles |17fl and laser-polarized 
Xe gas using NMR techniques |lS|| . 

In many of the nonequilibrium systems discussed 
above, the underlying stochastic process ip(t) is nonsta- 
tionary. For example, the two time correlation function 
C(t\,t 2 ) = (ip(ti )ip (£2)) for the diffusion equation de- 
pends on the ratio of the two times t\ and t 2 , and not 
on their difference p|. Persistence in such systems typ- 
ically decays as a power law P(t) ~ t~ 9 at late times 
t. The exponent 9, called the persistence exponent, is 
believed to be a new exponent and is apparently un- 
related to the usual dynamical exponents that charac- 
terize the decay of n-point correlation functions with fi- 
nite n. Persistence has also been studied for station- 
ary processes [jl9|,^0| such as a stationary Gaussian pro- 
cess characterized by its two time correlation C(t\,t 2 ) 
which is invariant under arbitrary time translation, i.e., 
C(ti + t , t 2 + t n ) = C(ti, t 2 ) for all t . In the stationary 
case, persistence between times t\ and t 2 typically de- 
cays exponentially, P(ti,t 2 ) ~ exp[— 9 s \t 2 — t\\] for large 
time difference [^0| . For some processes such as the diffu- 
sion equation the nonstationary problem can be mapped 



onto a corresponding stationary one || and the exponent 
of the nonstationary process becomes identical to the 
inverse decay rate 9 S of the corresponding stationary pro- 
cess Q. Despite many theoretical studies of either or 
6 S , exact results are known only in relatively few cases 
pl[ . The basic difficulty in computing either of them can 
be traced back to the fact that the underlying stochastic 
processes in both cases are usually non-Markovian 

In this paper we study the persistence in stochastic 
processes that are stationary under translations in time 
only by an integer multiple of a basic period (without 
loss of generality, this period may be chosen to be 1). 
Throughout this paper we will refer to such processes 
as SIS (stationary under integer shifts). For example, 
a Gaussian stochastic process will have the SIS prop- 
erty if its two-time correlation function C(ti, t 2 ) satisfies 
C(ti + n, t 2 +n) = C{t\,t 2 ) for all integer n. Such pro- 
cesses appear in many physical situations. For example, 
in weather records, there is an underlying non-random 
periodic forcing ( the motion of earth round the sun), 
which makes the stochastic process not truly stationary 
in time. In nonlinear systems, even if one can filter out 
the periodic component, the properties of the filtered sig- 
nal (say variance) would still be expected to show a pe- 
riodic variation with time. It seems worthwhile to study 
in more detail persistence in such SIS processes. 

When one wants to study the persistence of SIS pro- 
cesses, the following question arises naturally: Is the 
probability P{t) that the process remains positive over 
the interval [0, t] is same as the probability P n that the 
process is positive only at all the n intermediate integer 
times between and tl In other words, is the persistence 
of a 'continuous' SIS process is the same as the persis- 
tence of the corresponding 'discrete' sequence obtained 
by measuring the process only at integer times? 

The question regarding the difference between 
'continuous-time' and 'discrete-time' persistence was first 
raised in Ref. for strictly stationary Gaussian pro- 
cesses, motivated from the observation that in experi- 
ment |H| as well as numerical simulation |g3[ of persis- 
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tence in diffusion equation, the measurements (whether 
the process is positive) are done only at discrete points 
(separated by a fixed time window of size A) even though 
the actual process is continuous. In Ref. |^2| it was shown 
that for general stationary Gaussian processes, the con- 
tinuous persistence decays as P(t) ~ exp[— 9 c t] for large 
t, where as the corresponding 'discrete-time' persistence 
(obtained from measuring the data only at the inter- 
mediate time points separated by a fixed A) decays as 
P n ~ exp[— Odin] where T — nA. In general, one would 
expect that the exponent 9 C is strictly greater than 9d for 
such a process since the process can change sign between 
two successive integer points. The exponent 9d was com- 
puted analytically in Ref. for a Gaussian stationary 
Markov process and was shown to depend continuously 
on the window size A. 

In this paper we study the 'continuous-time' versus 
'discrete-time' persistence for SIS processes. We restrict 
ourselves to the study of discrete-time persistence only 
when the measurement points in time are integers. This 
is a natural choice since the process is stationary only 
under integer shifts in time. We show that unlike strictly 
stationary processes, in this case, the two exponents 9 C 
and 9d can be equal under certain conditions. This 
discrete-time persistence P n of a sequence is also rather 
interesting from a purely mathematical point of view, es- 
pecially when the underlying process tp(t) is Gaussian. 
In that case, calculating P n becomes the problem of cal- 
culating the probability that a set of n Gaussian random 
variables with a specified joint probability distribution 
are all positive. This 'one-sided barrier' problem has re- 
mained popular in the applied mathematics literature for 
many decades jl9|,^0],^4| . But the number of cases where 
this probability can be explicitly calculated for large n 
remains rather small p4j . 

In this paper, we study a specially simple case of 
a continuous-time stochastic process tp{t) which is ob- 
tained as a local smearing of a sequence of independent 
identically distributed random variables via a smearing 
function f(t). This process becomes, by construction, a 
stochastic process whose probability distribution is in- 
variant under time translations by integers, i.e. a SIS 
process. We construct examples where 9d < 9 C , and also 
construct a family of smearing functions f(t) for which 
Od = 9 C . We provide a physical example namely the diffu- 
sion equation on a hierarchical lattice where the diffusion 
field is a Gaussian stochastic process with the SIS prop- 
erty and we compute the corresponding smearing func- 
tion exactly. We then determine exactly the exponent 9d 
for a specific case when the correlations in the discrete se- 
quence are nonzero only for consecutive values. We find 
that in this case, the exponent 9d depends continuously 
on the value of the correlation. 

The paper is organized as follows. In section II, we 
give some examples of continuous SIS processes where 
9 C is strictly greater than 9d, and also construct a class 



of processes for which 9 C = 9d- In section III, we pro- 
vide a physical example namely diffusion equation on 
a hierarchical lattice where the diffusion process shows 
log-periodic oscillations. After rescaling, and a change of 
variables from time t to log(t), we get a stochastic process 
that has the SIS property. In section-Ill, we introduce 
a special sequence for which the persistence exponent 9d 
can be computed exactly. Section IV contains a summary 
of our results. 



II. CONTINUOUS-TIME VERSUS 
DISCRETE-TIME PERSISTENCE 

In this section we discuss the conditions under which 
the continuous-time persistence of a SIS process is the 
same as the discrete-time persistence of a correspond- 
ing sequence obtained by measuring the process only at 
the integer points. As noted earlier, in general, we expect 
that 9 C > 9d- Consider a stochastic process ip(t), which is 
known to be positive at the integer points t = 1, 2, . . .N. 
Now consider the conditional distribution of ip(ti) at 
some non-integer point t\ lying within the interval [1, N]. 
This conditional distribution is a gaussian whose width 
is independent of the values of ip(t) at the known inte- 
ger points. If this variance is finite, even in the limit 
when TV is large, clearly, there will be finite probability 
of i/j(ti) becoming negative. Thus one would expect to 
get 9 C = 9d only if the conditional variance of tpi^i) tends 
to zero for large N. In such a process, one should be able 
to determine ip(t) for all real t, if one knows its value at 
all integer points. 

This suggests the following construction of ip(t) : we 
consider a sequence of independent random variables 
{4>(n)} having zero mean, where n G { — oo,+oo}, and 
define a stochastic process ip(t) by the convolution 

m= E (i) 

n— — oo 

Knowing ip(t) at all integer points, one can expect to 
determine uniquely the constants {4>(n)} by solving cou- 
pled linear equations, which then determines ip(t) for all 
real values of t. 

The behavior of this process depends only on the 
smearing function f(t). In the following, we shall assume 
that f(t) has some good properties, i.e. is a non-negative 
unimodal function of t, which decreases sufficiently fast 
for large \t\. By a shift of the origin of time t, and rescal- 
ing tp(t), we can assume that the maximum of f(t) occurs 
at t = and /(0) = 1. 

What is the class of functions f(t) such that 9 C equals 
9dJ This class is not easy to characterize directly. A 
simple example illustrates this point clearly. Consider 
the simple case of triangular function 
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f(t) = 1 - \t\/a, for \t\ < a (2) 
= 0, otherwise. (3) 

In this case, ip(t) is a piece-wise linear function of t. 
If a < 1/2, we have intervals in which 4>(t) is identi- 
cally zero. If, however, we define persistence probabil- 
ity as the probability that the function does not change 
sign upto time £, it is clear that for all a < 1, we have 
9 c = 9 d = log2. 

We now show that 9 C ^ 0j if a > 1. For this purpose, 
it is sufficient to show that there are sequences {4>n} such 
that the corresponding ^-process is positive at all inte- 
ger points, but takes negative value for non-integer t. As 
such events would occur with non-zero probability along 
the sequence, 9 C > 64. 

Let n be the integer just below a. We consider a peri- 
odic sequence of </>„ with 

cf>i = 1, if i = mod(2n + 1); 

= — c, if i = n or n + 1 mod (2n + 1); 

= otherwise. (4) 

Then it is easy to see that if we choose c such that 
(a - n) > (2a - l)c > (a — n — 1/2), then ip(t = i) is 
positive for all integers i, but ip(t = n + 1/2) is nega- 
tive. Clearly, these signs are not changed if all </>'s devi- 
ate from these values by sufficiently small amount. Then 
such sequences ( of finite length) will occur with non-zero 
frequency, and hence for any a > 1, 9 C is strictly greater 
than 9d- 

However, there are functions /(£) for which 9 C = 64. 
The simplest example of this class is f(t) = exp(— |£|/a). 
In this case, it is easy to see from Eq. (|l|) that i])(t) at 
any non-integer point t can be expressed as a positive 
linear combination of its value at the two nearest integer 
points, so that for all £ = n + St with < St < 1 we have 

ip(n + St) = [sinh — tb(n + 1) (5) 
a 

+ sinh ^ ~ 5t ' Mn)]l [sinh(l/a)1 (6) 
a 

Thus, if ip(n) and ip(n + 1) are positive, Eq. (Q) implies 
that tp(t) is positive for all n < t < n + 1. Hence one 
gets, 9 C = 9 d - 

The example above can be generalized. For example, 
one can introduce a two parameter family of functions, 
f(t) = exp(ftii) for t < 0, and /(t) = exp(— re 2 i) for t > 
with Ki > and K2 > and are not equal in general. In 
fact, one can even introduce two arbitrary periodic func- 
tions g\ (t) and gi (t) ( with period 1), and take 

f(t) = exp( Kl t-g 1 (t)), for t < 0, (7) 
= exp(-/c 2 i - ffa(*))) f°r ^>0 (8) 

without destroying the equality of 9 C and One 
only has to impose some conditions on g\(t) and 52 (i) 



to ensure that /(£) is unimodal. Effectively, we can 
take any unimodal function f(t) defined in the interval 
— 1 < t < 1, and extend it to all real t using the conditions 
fit - 1) = e" Kl /(t), for i < 0, and f(t + 1) - e^/W 
for t > 0, to get a function /(i) for which 9d and # c are 
equal. 

III. PERSISTENCE IN DIFFUSION EQUATION 
ON A HIERARCHICAL LATTICE 

A simple example of a physics problem where functions 
of the type given by Eq. (^) show up is the persistence of 
a diffusion field on a hierarchical lattice. The lattice may 
be thought of as a line having N — 2 n sites labelled by 
an n-bit binary integer i,0 < i < N — 1 p5|p6[ . Wc 
define the ultrametric distance between two sites i and 
j as d, if the binary integers denoting i and j differ at 
the n — d + 1 bit counting from the left. Thus we have 
d = 1 between sites 2 and 3, but d = 4 between sites 7 
and 8. , At each site i, we have a real variable VK*)- At 
time t = 0, the fields at different sites are assumed to 
be independent identically distributed random variables 
(say gaussians of mean zero, and variance 1). The fields 
are assumed to evolve in time by the deterministic 
equation 

N-l 

j-m = Y, K ^3)-m] (9) 

Here the spring constants Kij are assumed to be func- 
tions of the distance dij between the two points. In the 
following, we shall assume that Kij = a~ diJ , where a is 
a constant > 1. 

The integration of the equations of evolution is made 
particularly simple by the hierarchical nature of the 
spring couplings. It is easily verified that we have 2™~ r 
independent eigenmodes of relaxation rate (a— l) _1 a _r+1 
(r = 1,2, . . . N - 1) satisfying 

j t Sf ] = (a- l)- 1 a-'"+ 1 sj r3 (10) 

where 

2 r - 1 -l 

4 )= E 0(j2 r - 2 1 - 1 - fc) - ^ - fc) (11) 

fc=0 

where j = 1 to 2™~ r . 

Expanding any particular say for i = 1, in terms 
of these eigenvectors, and we get, 

n-l 

Mt) - E 2 ~ r/2ex P ["(« - l)- 1 o- r+1 i] #r), (12) 

r=l 

where </>(r)'s are i.i.d. Gaussian variables of zero mean 
and unit variance that characterize the initial condition. 
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This formula for the hierarchical model may be compared 
with the corresponding formula one writes in the Eu- 
clidean space in d— dimensions 

p oo 

t/j(r = 0)= dkcxp{-k 2 t)ri(k) (13) 
Jo 

where rj(k) are white-noise process with variance 

< rj(k)i](k') >= 5 ktk ,k d - 1 (14) 

We eliminate the time variable t in terms a logarith- 
mic time variable r using the identification a T = i, and 
change ip(t) by a change of scale, tp(r) = [a T ' 2 ]ip(t — 
(a — l)a T ). Then we have 

oo 

V>(r) = f M ex P ha"] a (T_r)/2 - (15) 

For large r, the summation over r can be extended from 
—oo to +oo, and the process ip(r) then becomes a Gaus- 
sian process with the SIS property, i.e., is stationary 
only under integer shifts in time and is obtained by local 
smearing of the discrete white noises </>(r)'s, 

E f(r-r)<f>(r), (16) 

r— — oo 

where the convolution function / (r) clearly goes to zero 
when r tends to ±oo. Thus, the problem of calculating 
the persistence exponent reduces that of calculating the 
exponents 9 C and 9d for a process defined by given con- 
volution function f(t) — exp(— a*)a'/ 2 . The origin of the 
SIS property here comes from the discrete scale invari- 
ance of the original model, which gives rise to log-periodic 
oscillations in the relaxation processes (27] . 

We have not been able to determine whether this func- 
tion fit), the exponents 9 C and 9d coincide, or are differ- 
ent. However, in a simple Monte Carlo realization of a 
sequence of 10 5 Gaussian variables {<fii}'s of zero mean 
and unit variance, we did not find any instance where 
the function ip(t) changed sign twice between two con- 
secutive integers. This indicates that these exponents, if 
not equal, are likely to be quite close to each other. 

IV. EXACT RESULTS FOR A SPECIAL CASE 

For a smearing function f(t) for which 9 C = 9d, the 
computation of the persistence exponent simplifies con- 
siderably, and reduces to its determination for a discrete 
sequence rather than a continuous process. But even 
then, the exponent 9d is quite nontrivial for an arbi- 
trary smearing function f(t). For calculating 9d, only 
the values of f(t) at integer points are relevant. In the 
following, we shall consider in detail the calculation of 
9d when only /(O) and /(— 1) are non-zero. This can be 



thought of a crude approximation to the smearing func- 
tion f(t) = cxp(— a')a*/ 2 , as in the diffusion equation 
on a hierarchical lattice, which decreases superexponen- 
tially for t > and only exponentially for t < for a > 1 . 
We will show below that the exact computation of 9d is 
nontrivial even for this toy smearing function since the 
resulting sequence is non-Markovian. 
In this special case, Eq. ( |l6| ) becomes 

ipi = 4>i +e<f> i+ i, i=l,2,...,n (17) 

where we shall assume that {<fii} are independent identi- 
cally distributed random variables, not necessarily Gaus- 
sian, each drawn from the same distribution p{(p). Here 
e is a mixing parameter. For convenience, we relabel the 
0's without any loss of generality to consider the follow- 
ing sequence, 

ipi = 4>i + € <t>i-ii i =1,2,. . . ,n. (18) 

For simplicity, we will assume that p{<p) is symmetric 
about the origin. The mean value of <fi is then zero. We 
now ask: what is the probability P„(e) that ipi, ip2, ■ ■ ■ ipn 
are all positive for a given e? 

We note that the variables ^'s are now correlated. 
The two point correlation function, C^.j = {i^itpj) can be 
easily computed from Eq. (|l~8|), 

C\ 3 = a 2 [(1 + e 2 )5 itj + ,u>, l ; + Sij-x)] , (19) 

where Si j is the Kronecker delta function and a 2 — 
/_ 4> 2 p{4>)d4>- Thus the parameter e serves as a measure 
of the correlation and it is this correlation that makes the 
calculation of P n (e) nontrivial for nonzero e. 

The sequence {i/j n } defined by Eq. (|l^) is non- 
Markovian in the sense that if only {ipn} are observed, 
and not the ^> n 's, ipn depends not just on the previous 
member of the sequence VVi-i, but rather on the whole 
history of the sequence. For example, from Eq. ( |l8| ) one 
can express ipn as, 

n-l 

^ n = ^{-^) k ~ 1 ^n-k + 4> n -e n <h (20) 
fc=0 

which demonstrates explicitly the history dependence of 
the sequence. For non-Markovian sequences, it is gen- 
erally hard to compute the persistence exponent. For- 
tunately progress can be made for this special sequence 
even though it is non-Markovian. 

In order to calculate P n (e), it is first useful to define 
the following probabilities, 

Qi(x) = / d^ o p(0o), 

J x 

pQC pQO pOC 

Q„(x)= / d(f> p((j) ) / <20i/?(0i) / d(p 2 p{(p2) ■ ■ ■ 

J x J — e0o J —e4>\ 




d^„_ip(^_l), n > 2. (21) 
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Using the definitions in Eq. (|l8|), it is then easy to see 
that the persistence P n (e) = Q n +i(—oo). This is due to 
the fact that for all the i/j/s in Eq. ([lj) to be positive, 
while (j)Q is free to take any value, <pi must be bigger than 
—e(f>o, (f>2 must be bigger than — e<j>i and so on. Differen- 
tiating Eq. ( pi] ) with respect to x, we get the recursion 
relation 



dQnjx) 
dx 



-p(x)Q n -x(-ex), n>l, (22) 



with Qo(x) = 1 and the boundary condition, Q n {oo) = 
for all n > 1. Let us define the generating function 



F(x,z) = J2Qn(x)z n . 



(23) 



^From Eq. (p2|), it follows that F(x,z) satisfies a first 
order non-local differential equation, 



dF(x, z) 
dx 



= -p{x)z [1 + F(-ex, z)] , 



(24) 



with the boundary condition, F(oo, z) = for any z. 
Once we know the function F(x, z), P n (e) can be ob- 
tained by evaluating the Cauchy integral, 



P n {e) = Q„+i(-oo) = — / 

2 ™ J Go 



1 f F( -oo,z) J 

^—dz, (25) 



over a contour Co encircling the origin in the complex z 
plane. 

Before proceeding to solve Eq. ( p3| ) , we make the sim- 
ple observation that, 



P n (e)=P n 



(26) 



true for any e. To see this, we first rescale the ipi vari- 
ables, ipi = ipi/e. Clearly the persistence of ip/'s is the 
same as that of the tpi 's. Dividing Eq. ( p"8"| ) by e, we 
see that in order for the ip/'s to be positive, we need to 
satisfy the conditions: (f>o > —(f>i/e, 4>\ > —fo/e, ■ ■ ■, 
4> n -\ > — ^n/e where <j) n can be arbitrary. Eq. (pq) then 



follows once we relabel i 



for all < i < n. Thus 



it is sufficient to compute -Pn(e) for e only in the range, 

— 1 < e < 1. Once we know this, P n (e) for |e| > 1 can be 
obtained from Eq. (^6|). 

Let us summarize our main results. We show that for 

— 1 < e < 1, fn(e) ~ exp(— 9(e)n) for large n, where 
6(e) depends continuously on e and also depends on the 
distribution p((f>). In contrast, at e = 1, the exponent 
6(1) = log^) is independent of the distribution p(<fi). 
The exponent 6(e) diverges as e — > — 1, indicating a faster 
than exponential decay of P n for large n. We show that 
P n (— 1) = l/(n+ 1)! exactly for all n > 1, again inde- 
pendent of the distribution p(4>). 



A. The case when e = — 1 

Let us first consider the case e = — 1. In this case, the 
Eq. (|24|) becomes local and can be easily solved by inte- 
gration. For symmetric p((f>) with zero mean, the exact 
solution is given by 



F(x,z) 



- 1 + exp 



1 



p(x')dx' 



(27) 



which satisfies the boundary condition F(oo, z) = 
for all z. Expanding the exponential in Eq. (|27| ) in 
powers of z and using the definition in Eq. (|23|), we 
find, Qn(x) = (5 — Jo p(x')dx') 71 /nl. Using the rela- 
tion P n = Q n+ i(—oo) and the normalization condition 
I- 00 P( x ')dx' = 1, we get 



Pn(-l) 



1 



(n + 1)! ' 



(28) 



for all n > 1. Remarkably P n (—1) is independent of the 
distribution p((f>) for all n > 0. 



B. The case when e = 1 

Next we consider the case e = 1. We first make a 
change of variable, u(x) = J p(cj))d<p. Let F(x, z) — 
F(u,z). Since p((f>) is symmetric about zero, u(— x) = 
—u(x) and hence F(—x, z) — F(—u, z). Using this in Eq. 
(p4|) with e = 1, we find 



dF(u, z) 
du 



l+F(-u,z) 



(29) 



where u varies from —1/2 to 1/2 and the boundary con- 
dition is, F(l/2, z) = Q for all z. Differentiating Eq. (§|) 
with respect to u we get a local second order differential 
equation 



d 2 F(u,z) 
du 2 



1 + F(u,z] 



(30) 



whose general solution is given by 

F(u, z) = -1 + \A(z) cos(zu) + B(z) sin(zu)] . (31) 

If this solution also has to satisfy Eq. (|2^), we have 
additionally B(z) = —A(z). The boundary condition 
_F(l/2, z) — determines A(z) and we finally get 



F(x,z) = -1 



cos (u(x)z) — sin (u(x)z) 
cos(z/2) — sin(2;/2) 



(32) 



Thus, F(-oo,z) = 2/[cot(z/2) - 1]. This function has 
poles at z = 7r/2 + 2mn, where m is an integer. One can 
then easily evaluate the contour integration in Eq. ( |25| ) 
and we get the exact expression, 
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2m7r) 



n+2 ! 



(33) 



valid for any n > 0. For example, by summing the series 
in Eq. © we find, P (l) = 1, P x (l) - 1/2, P 2 (l) = 1/3, 
^3(1) = 5/24, etc. which can also be verified by per- 
forming the direct integration in Eq. (pl|). The remark- 
able fact is -Pn(l) is universal for all n > in the sense 
that it is independent of the distribution p(4>), as in the 
e = — 1 case. Clearly the leading asymptotic behavior is 
governed by the m = term in Eq. ( j33| ) and we get, 
F n (l) - exp(-9n) for large n, with 0(1) = log(-7r/2). 
Clearly the exponent 9(1) is also universal. 

Interestingly, P n (l) is related to the fraction of 
metastable states in an Ising spin glass on a 1- 
dimensional lattice of n sites at zero temperature [ p8| . 
Consider the spin glass Hamiltonian on a chain, H = 
— Ji,i+i s i s i+i where Si — ±1 are Ising variables and 
the bonds J^+i's are independent and identically dis- 
tributed variables each drawn from the same symmetric 
distribution with zero mean. Out of the 2" number of to- 
tal configurations, how many are metastable with respect 
to single spin flip Glauber dynamics at zero temperature? 
A configuration is metastable at zero temperature if the 
energy change AE, t = 2s l [J i - M s i _i + J ii+1 s i+ i] > 
due to the flip of every spin Si . Defining the new variable 
4>i = 2 Ji,i+iSiSi+i, we see that the variables 4>iS are also 
independent and identically distributed and the proba- 
bility that a configuration is metastable is precisely the 
probability that the variables, ipi = 4n + 4>i-i are posi- 
tive for each i. This is precisely P n [X) as computed in the 
previous paragraph. We note that the average number of 
metastable configurations (N s ) for the 1-d spin glass was 
computed exactly by Derrida and Gardner |2{| by a dif- 
ferent method and they found (N s ) ~ (4/7r) ra for large n. 
Thus the fraction of metastable configurations scales as 
{N s )/2 n ~ (ir/2)~ n , in agreement with our exact result 
forP„(l). 



C. The case — 1 < e < 1 

We now turn to the range, — 1 < e < 1. In this range, 
we were unable to calculate P n (e) exactly for arbitrary 
distribution p(4>). However progress can be made for the 
uniform distribution, 



p(<f>) 



for 



1<6<1 



— 0, otherwise. 



(34) 



For this case, it follows from Eq. ( p4| ) that F(x, z) is 
independent of x for x < — 1 and hence, F(— 00, z) = 
F(-l,z). Similarly, F(x,z) = for all x > 1. In the 
range, — 1 < x < 1, we expand F(x,z) = b m (z)x m 
in a power series in a;. Subsituting this series in Eq. (24), 



we get the recursion relation, b m = — b m -\z(— e) m ~ 1 /2m 
for all m > 1. Thus the function F(x, z) can be expressed 
completely in terms of only b§(z) which is then deter- 
mined from the boundary condition, F(l,z) = 0. This 
determines F(x, z) completely in the range —1 < x < 1 



and we find, F(x, z) = —1 



where 



/(*) = E 



(_T)Mrn+lV2 _ 1)/2 



m=0 



(35) 



Using F(— 00, z) = F(—l, z), we finally get 

/(-*) 



F(-oo,z) = -1 



(36) 



where f(z) is given by Eq. (|35|). We note that the series 
in Eq. ( pq ) and hence in Eq. ( |36| ) is convergent for all 
z as long as — 1 < e < 1. In fact, for e = 1, it is easy to 
see that Eq. (||) gives F(-oo,z) = 2/[cot(z/2) - 1] as 
before. 

Substituting Eq. (|6|) in the expression of P n (e) in Eq. 
( p5| ) , we find that the leading asymptotic decay of P„ for 
large n is governed by the pole of F(— 00, z) that is closest 
to the origin. From Eq. (|36"|), the poles of F(— 00, z) in 
the z plane are precisely the zeroes of the function f(z) 
in Eq. ( |35| ) in the z plane. In particular, P„(e) ~ z + ~ n 
for large n where z+ is the zero of f(z) in Eq. (|3j]) closest 
to the origin. The persistence exponent is then simply, 
9 = log(z + ). Let us first consider a few special cases. For 
e = 1, we find from Eq. (|35|), f(z) — cos(z/2) - sin(z/2) 
indicating z + = 7r/2, a result we already obtained. For 
e = 0, we find from Eq. (p5|), f(z) = 1 — z/2, indicating 
z + = 2, as expected for the persistence of uncorrelated 
variables. As e — > — 1 + , the function j(z) is Eq. j3^ ) ap- 
proaches to, f(z) —> exp(— z) indicating z + —> 00. Indeed 
putting e = — 1 + S in Eq. (|35|), it is easy to see that, 
z + « \/$>/S as S — > 0. Thus the persistence exponent 
diverges as, 9 m log[-\/8/ (1 + e)] as e — > — 1. 

For other values of e in the range — 1 < e < 1, it is easy 
to evaluate z+ to any arbitrary accuracy from Eq. ( |35| ) 
using Mathematica. The exponent 9 = log(z+) for some 
representative values of e in the range — 1 < e < 1 are 
listed in Table 1. The exponent 9(e) increases monotoni- 
cally as e decreases from +1 to —1, diverging as e — > — 1. 
For |e| > 1, the exponent is determined from the relation, 
9(e) = 9(1 je). Thus in the whole range, —00 < e < 00, 
the exponent 9(e) is a nonmonotonic function of e. As 
e varies from — 00 to 00, 9(e) increases monotonically 
in the range [—00, —1], then decreases monotonically in 
[—1, 1] followed by a further monotonic increase in the 
range [l,oo]. The slowest decay of P n occurs at e = 1, 
where 9(e) is minimum and given by the universal value, 
0(1) = log(7r/2). 



G 



e 


6> 


1.0 


loe:f7r/2 > ) = 4515 


3/4 


0.4690 . . . 


1/2 


0.5155 


1/4 


0.5882... 





log(2) = 0.6931 . . . 


-1/4 


0.8465 . . . 


-1/2 


1.0906... 


-3/4 


1.5686... 



Table 1. The exponent 9(e) for some representative val- 
ues of e in the range, — 1 < e < 1 in the case of the 
uniform distribution in Eq. (|3^). 

Except at e = 0, 1 and —1, the exponent 9(e) is nonuni- 
versal in the sense that its value depends on the details 
of the distribution p(4>). To see this explicitly, we now 
compute 0(e) perturbatively for small e. We expand the 
right hand side of Eq. ( p4| ) upto order e and then solve 
the resulting local differential equation exactly to deter- 
mine F(x,z) upto O(e). Taking x — > — oo limit in the 
expression of F(x, z), we find 



F(-oo,z) 



2z 



(37) 



where c = J 
to the origin is given by. 



[2 - z - 2cp(0)ez 2 ] ' 

From Eq. (|37|), the pole closest 

(38) 



z + = 2[l-4cp(0)e + O(e 2 )] . 



From Eq. (|25j), it then follows that P n (e) ~ z+~ n for 
large n. Hence 9(e) = log(z + ) = log(2) - 4cp(0)e + 0(e 2 ) 
and is clearly nonuniversal, as seen from the nonuniver- 
sality of the O(e) term in the previous equation. For 
example, for the uniform distribution in Eq. (|34|) , we get 
9(e) = log(2) - e/2 + 0(e 2 ). On the other hand for the 
Gaussian distribution, p(4>) = (2n)^ 1 / 2 exp(— <j) 2 /2), we 
get 9(e) = log(2) - 2e/n + 0(e 2 ). 



V. CONCLUSION 

In summary, we have discussed persistence in stochas- 
tic processes that are stationary only integer translations 
of time. Such a process can be explicitly constructed by 
smearing independent noises with a convolution function 
f(t). A physical example of such a process is provided by 
the diffusion field on a hierarchical lattice for which we 
have computed the smearing function f(t) exactly. How- 
ever, we could not compute the persistence exponents 9 C 
or 9d in this case. We showed that under certain condi- 
tions, the continuous-time persistence of such a process 
reduces to the persistence of a discrete sequence obtained 
by measuring the process only at integer times. We have 
constructed a specific non-Markovian sequence where the 
smearing function is nonzero only at two consecutive in- 
teger points leading to nonzero correlations only between 



consecutive values of the sequence and computed the per- 
sistence exponent 84 exactly for this sequence. The ex- 
ponent 9d depends continuously on the strength e of the 
correlation. Remarkably for e — 1 and e — — 1, the per- 
sistence becomes universal. For e = 1, we have shown 
an interesting connection between the persistence of this 
sequence to the average fraction of metastable states in 
a one dimensional spin glass. 

The class of functions f(t) for which we could show 
that 9 C = 9d is perhaps not the most general. A pre- 
cise characterization of this class seems like an interest- 
ing problem. Calculation of the persistence for SIS pro- 
cesses, or sequences, with correlations extending beyond 
nearest neighbors may be possible in some special cases, 
and would help understand the general question about 
the dependence of the persistence exponent on the cor- 
relations in the sequence. 

We thank A. Lefevre, D.S. Dean, V. S. Borkar and A.J. 
Bray for useful discussions, and M. barma for a critical 
reading of the manuscript. 
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